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<^ : Abstract 



Two-body charmless hadronic B decays involving a light tensor meson in the final states are 
studied in the perturbative QCD approach based on kx factorization. From our calculations, 
we find that the decay branching ratios for color allowed tree-dominated decays B — t- a2'7r^ and 



,^ ! B — )■ On vr+ modes are of order 10 ^ and 10 ^, respectively. While other color suppressed tree- 

f~^ ' dominated decays have very small branching ratios. In general, the branching ratios of most 

D ■ _^ _s 

r~| ' decays are in the range of 10 to 10 , which are bigger by one or two orders of magnitude than 



those predictions obtained in Isgur-Scora-Grinstein-Wise II model and in the covariant light-front 
approach, but consistent with the recent experimental measurements and the QCD factorization 
calculations. Since the decays with a tensor meson emitted from vacuum are prohibited in naive 
factorization, the contributions of nonfactorizable and annihilation diagrams are very important to 
these decays, which are calculable in our perturbative QCD approach. We also give predictions to 
the direct CP asymmetries, some of which are large enough for the future experiments to measure. 
k><( ■ Because we considered the mixing between /2 and /2, the decay rates are enhanced significantly 

C^ . for some decays involving /2 meson, even with a small mixing angle. 
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I. INTRODUCTION 



In the quark model, all kinds of mesons are classified by the spin-parity quantum numbers 
J^ . For example, J^ = 0~ denotes pseudoscalar mesons and J^ = 2"*" represents tensor 
mesons. The p-wave tensor mesons that we study in this paper include isovector mesons 
02(1320), isodoublet states K;{U30) and two isosinglet mesons /2(1270), /2(1525) |l|, bj. 
For these nine tensor mesons, both orbital angular momentum and the total spin of quarks 
are equal to 1. Because of the requirement of the Bose statistics of the tensor meson, the 
light-cone distribution amplitudes of tensor mesons are antisymmetric under the interchange 
of momentum fractions of the quark and anti-quark in the flavor SU(3) limit (3|, |j]. 

Recently, several experimental measurements about charmless B decay modes involving a 
light tensor meson (T) in the final states have been obtained J5Nl8|. These decays have been 



studied in the naive factorization approach 19l-|27l|. with which it can be easily shown that 



\i 



22 



23| . The factorizable 



(0 I j'' I T) = 0, where jf" is the {V ± A) or {S ± P) current 
amplitude with a tensor meson emitted vanishes. So these decays are prohibited in the naive 
factorization approach. The branching rations predicted in the naive factorization approach 
are too small compared with the experimental results, which implies the importance of 



nonfactorizable anc 
approach analysis 



annihilation type contributions. The recent QCD factorization (QCDF) 
41 proved this. It is worth of mentioning that the perturbative QCD 
(PQCD) approach [28|, |29| is almost the only method to calculate these kinds of diagrams, 
without fitting the experiments. 

In this work we shall study charmless -B„(d) -^ PT decays in the perturbative QCD 
approach based on the k-r factorization. Due to the heavy mass of B meson, the two light 
mesons decayed from the B meson are moving very fast in the rest frame of B meson. The 
light quarks in the final state mesons are all coUinear; while the light spectator quark from 
B meson is soft. Therefore there must be a hard gluon to kick the light spectator quark in 
the B meson to form a fast moving light meson. In this case, the hard process dominates 
the decay amplitude, which make it perturbatively calculable. By keeping the transverse 
momentum of quarks, the end point singularity in the coUinear factorization can be elim- 
inated. Double logarithm appears in the QCD radiative corrections due to the additional 
energy scale introduced by the transverse momentum. By using the renormalization group 
equation, the double logarithm can be resumed and leads to the Sudakov factor, which ef- 



fectively suppresses the endpoint contribution of the distribution amphtude of mesons in 
the small momentum region to make the perturbative calculation reliable. The annihilation 
diagrams can also be perturbatively calculated in the PQCD approach, which is proved to be 



the dominant strong phase in B decays for the direct CP asymmetry SOj. Phenomeno logi- 



cally, the PQCD approach has successfully predicted the direct CP asymmetry in hadronic 



B decays [30| and the branching ratios of pure annihilation type B decays J31|. 

This paper is organized as follows. In Sec. II, we present the formalism and wave func- 
tions of the considered B meson decays. Then we perform the perturbative calculations for 
considered decay channels with the PQCD approach in Sec. III. The numerical results and 
phenomenological analysis are given in Sec. IV. Sec.V contains the main conclusions and a 
short summary. Finally, Appendix A contains input parameters and distribution ampli- 
tudes used in this paper and Appendix B gives various functions that enter the factorization 
formulae in the PQCD approach. 

II. FORMALISM AND WAVE FUNCTIONS 



The related weak effective Hamiltonian -ffe// 32| for charmless b — > d{s) transitions can 
be written as 

r ( ^ 10 >| 

Heff = -^\H a(^)K*6KDor(/i) - v,iVtD Yl cAiAOM , (1) 

I j=i i=3 j 

where Vub-, Vud, Vtb and Vto are CKM matrix elements, D denotes the light down quark d or 
s, and Ci(^j){fi) are Wilson coefficients at the renormalization scale /i. Oi(^j){fi) are the well 



known effective tree (penguin) operators 32|. 



The non-leptonic B meson decays involve three energy scales, including the electroweak 
scale M\Y, b quark mass scale Mb and the factorization scale v AM^, where A = Mb — mi,. 
When the energy scale is higher than the W boson mass M^/, the physics is the electroweak 
interaction which can be calculated perturbatively. The physics from Mw scale to Mb 
scale is described by the Wilson coefficients of effective four quark operators, which is the 
resummation of leading logarithm by renormalization equations. The physics between Mb 
scale and the factorization scale is calculated by the hard part calculation in the PQCD 
approach. The physics below the factorization scale is described by the hadronic wave 
functions of mesons, which are nonperturbative but universal for all decay processes. 



In the PQCD approach, the decay amphtude can be factorized into the convolution of the 
Wilson coefficients, the hard scattering kernel and the light-cone wave functions of mesons 
characterized by different scales, respectively. Then, for B — ?■ M2M3 decays, the decay 
amplitude is conceptually written as the convolution. 



A ^ / dxidx2dx3bidbib2db2b3db3 

xTr[C(t)$B(xi, bi)<^MA^2, b2)<^M,{^3, b^)H{x,, h, t)St{xi)e-''% (2) 

where Xi is the longitudinal momentum fractions of valence quarks, bi is the conjugate space 
coordinate of the transverse momentum kiT of the light quarks, and t is the largest scale 
in function H{xi,bi,t). By using the renormalization group equations, the large logarithms 
ln{my\r/t) are included in the Wilson coefficients C{t). By the threshold resummation, 
the large double logarithms In^ Xj are summed to give St{xi) which smears the end-point 
singularities on Xi 33|. The last term, e"'^^*-', is the Sudakov factor which suppresses the 



soft dynamics effectively [3J]. Thus it makes the perturbative calculation of the hard part 
H applicable at intermediate scale, i.e., rriB scale. 

We will work in the B meson rest frame and employ the light-cone coordinates for mo- 
mentum variables. So the B meson momentum is chosen as Pi = ^(1, 1, Ot)- For the 
non-leptonic charmless B — ?■ M2M3 decays, we assume that the M2(M3) meson moves in 
the plus(minus) z direction carrying the momentum P2{P3). Then the momenta are given 
by 

P2 = ^(1 - rl rl 0^), P3 = ^{rl 1 - r,^ 0^), (3) 



where r2 = -^^ and r^ = --^^, The (light-) quark momenta in B , M2 and M3 mesons are 



— - and rq = — ^ 
defined as fci, k2 and k^,, respectively. We choose 

ki = (XiP+, 0, kir), k2 = {X2P^, 0, k2T), ks = (0, 0:3^3", k3T). (4) 

For a tensor meson, the polarization tensor ef^iy{X) with helicity A can be constructed via 



Q 



the polarization vectors of a vector meson [3|, |J]. They are given by 

e^"'{±2) = e(±l)^e(±l)^ 

e^^{±l) ^ A/|[e(±ire(Or + e(0)M±i: 



e'^^(o) = ]ll[<+ir<-ir + e{-ire{+ir] + ^Jle{ore{or. (5) 



With the tensor meson moving on the plus direction of the z-axis, the polarization vectors 
of the vector meson are chosen as 



e'^(O) = -^^(ko + h, h - h, 0, 0), e'^(±l) = -^(0, 0, 1, ±t) 



(6) 



where fco denotes the energy and k^ is the magnitude of the tensor meson momentum in the 
B meson rest frame. The polarization tensor satisfies the relations 3|, |4| 



eA'-(A) = e-'^'iX), 



<(A) = 0, 



e^^(A)P^ = e'^'^(A)P, = 0, e^,(A)(e'^^(A'))* = ^av- 



(7) 



In the following calculation, we define a new polarization vector e^ for the considered tensor 
meson for convenience 



6t(A) = — e^,(A)P]^, 
niB 



(8) 



which satisfies 



e.,(±2) = 0, e.,(±l) = ^^^l^^^^i^, e.,(0) 



e(0) ■ PBe(O) 



(9) 



\/2mB ■ ' ■ rriB 

One can find that the new vector er is similar to the polarization vector e of a vector meson, 
regardless of the related constants j^ . 

In the PQCD approach, we should choose the proper wave functions for the B meson 
and light mesons to calculate the decay amplitude. Because the B meson is a pseudoscalar 
leavy meson, the two structure (7/^75) and 75 components remain as leading contributions 
2|. Thus the B meson wave function $5 is written as 



<l>B = -^[(f + ms)750B(a;)] 



(10) 



For the distribution amplitude, we can choose 

(f)B{x, h) = Nbx'^{1 — x)^exp 



1 /uibx 
2\^ 



colb' 



(11) 



where Nb is the normalization constant. 

For the light pseudoscalar meson (P), the wave function is generally defined as 

i 



$P X 



Vg 



75 {r<Pi{x) + ml^<p'p{x) + m^ifi^ - l)0?(x)} 



(12) 



where (^p' ' and ttiq are the distribution amphtudes and chiral scale parameter of the 
pseudoscalar mesons, respectively, x denotes the momentum fraction carried by the quark 
in the meson, and n = (1, 0, 0) and v = (0, 1, 0) are dimensionless light-like unit vectors 
pointing to the plus and minus directions, respectively. 

The wave functions for a generic tensor meson are defined by (2] 

1 



*^^7! 



F ■ V 



Here n is the moving direction of the tensor meson, and v is the opposite direction. We 
adopt the convention e°^^^ = 1. The vector e, = '''^ rriT is related to the polarization 



fir 

m 



tensor. The distribution amplitudes can be given by |2|-|4| 

^t{^) = 2^9?ix), Mx) = ^^^^9^±i^)- (14) 

The asymptotic twist-2 distribution amplitude is given by 

0ll^^(x) = 30x(l - x){2x - 1). (15) 

The twist-3 distribution amplitudes are also asymptotic and the forms are chosen as |2H4| 

/iW(x) = ^(2a;-l)(l-6a; + 6x2), h\'\x) = 15a;(l - a;)(2a; - 1), 

g^_^\x) = 20a;(l-x)(2x-l), g^l\x) = 5(2a; - 1)1 (16) 

III. PERTURBATIVE CALCULATION 

In this section, we will calculate the hard part H{t), which includes the effective four 
quark operators and the necessary hard gluon connecting the four quark operator with 



the spectator quark [36|. There are 8 types of diagrams contributing to the B — )■ PT 
decays, shown in Fig.l. From the first two diagrams of Fig.l, (la) and (lb), by perturbative 
QCD calculations, we gain the decay amplitudes for factorizable emission contribution. For 



T** <— 



(la) 



(16) 



(Ic) 



(Id) 





FIG. 1: Diagrams contributing to the B — )• PT decays, with a pseudoscalar meson emitted. 



(y — A)iy — A) current, the amplitude is written as. 



JG 

X {[0T(a;3)(a;3 + 1) - (^tI^^s) + 0T(x3))rT(2x3 - 1)] hef{xi,X2,MM)Ee}{ta) 

+ [2rT0T(a;3)] /ie/(a;3, Xi, 63, &l)^e/(tb)} , (17) 

where tt = ^, and Cp = |. /p is the decay constant of the pseudoscalar meson. The 
function hef, ta^b and E^j can be found in Appendix B. Form Eq{T71 we can obtain the 
{T\V — A\B) transition form factor in the PQCD approach. 

The operators O5, Oq, O7, and Og have the structure of (V — A){V + A). In some decay 
modes, some of these operators will contribute to the decay amplitude. Because only the 
axial part of {V + A) current will contribute to the pseudoscalar meson production, we have 






-A^i 



(18) 



In some cases, in order to get the right color structure, we must do a Fierz transformation 



for these operators. So we obtain [S — P){S + P) operators from {V — A)(y + A) ones. The 
decay amphtude is, 

AeT = ^Q\ h^CFfp'n"m'^B / dxidxs / bidbibsdbs ■ (f)B{xi,bi) 

X {[0^(3:3) +rTi<pTix3)ix3 + 2) - 0^(3:3) xg)] ro/ie/(a:i, X3, 61, 63)^6/ (ta) 

+ [2rTro0T(a^3)] hef{x3, Xi, 63, bi)Eef{tb)} , (19) 

where Tq = itLq /rriB- 

For the non-factorizable diagrams Fig.(lc) and (Id), the amphtudes involve all three wave 
functions. The integration of 63 can be performed through 6 function 6{bi — 63), leaving only 
integration of 61 and 62- For the (V-A)(V-A), (V-A)(V+A) and (S-P)(S+P) type operators, 
the amplitudes are 

32 f^ f°° 

M^T = — CpT^ra^B I dxidx2dx-i I bidbib2db2(j)Bixi,bi)(f)p{x2) 
3 Jo Jo 

X { [(Pt{x3){x2 - 1) + {(Pri^s) - <pTi^3))rTX3] 

■ henfiXl, 1 - X2, X3, 61, 62)^en/(ic) 

+ [(Pt{x3){x2 + X3) - (0^(3:3) + (l)T{x3))rTX3] 

■ henf{Xl,X2, X3, 61, 62)^en/(td)} , (20) 

/A^x = — — Ci?7rrom^ / dxidx2dx3 I bidbib2db2(j)Bixi,bi) 
3 Jo Jo 

X { [(PI{x2){(Pt{x3){x2 - 1) + rT{(pT{x3){-X2 + OJg + 1) + (Pt{^3){x2 + X3 - 1))) 
+<Pp{<Pt{x2 - 1) + rT{<l)Ti^3)ix2 - Xg - 1) - 0^(X3)(X2 + X3 - 1)))] 
■henf{Xi, 1 — X2, Xs, 61, b2)Egnf(tc) 

+ [(pp{x2){(pT{x3)x2 + rr(0T(a^3)(a;3 - 2:2) + 0T(^3)(a;2 + 3:3))) 

+(l)p{rT{<pTix3){x3 - X2) + <pTix3){x2 + X3)) - (I)t{x3)x2)] 

■ henf{Xi,X2, X3, 61, 62)-Ee„/(td)} , (21) 

00 /»1 /»oo 

A^e-r = ~";rC'F7rm^ / dxidx2dx3 / bidbib2db2(j)Bixi,bi)(j)p{x2) 
3 Jo io 

X { [0T(a;3)(a;2 - 3:3 - 1) + (0^(0:3) + 0r(x3))rTa;3] 

■ henf{Xi, 1 — X2, X3, 61, b2)Eenf(tc) 
+ [(Pt{x3)x2 + {4>T - 0T)^Ta;3] 

■ KnAxx,X2, X3, 61, &2)^en/(^d)} • (22) 



The factorizable annihilation diagrams Fig.(le) and (If), the three kinds of decay amph- 
tudes for these two diagrams are 

^aT = ^\j ^CpfBT^rn^ I dx2dx3 / b2db2hdb3 

X { [2(/)^{x2)rTro{(l)T{x3){x3 - 2) - (I)t{x3)x3) - (PJ>{X2)(PT{X3){X3 - 1)] 

■haf{x2, 1 - X3, 62, b3)Eaf(te) 

+ [2(j)'T{x3)rTro{(j)l{x2){x2 - 1) + (j)p{x2){x2 + 1)) - (j)i{x2)(j)T{x3)x2] 

■ haf{l - X3, X2, 63, b2)Eaf{tf)} , (23) 

A'J = -A^aT. (24) 

A^^ = IQ\ -CpfBT^m^B I dX2dX3 / 62(i62&3'^&3 

X { [20^(x2)0T(a;3)ro + 0p(x2)(0^(x3) + (t)T{x3))rT{x3 - 1)] 

• haf{x2, 1 - X3, 62, b3)Eaf{te) 

- [x2rQ(t)T{x3){(ty^p{x2) - (t)p{x2)) + 20^(a;2)0T(a^3)^T] 

• /la/(l - 3:3, X2, 63, b2)Eaf{tf)} . (25) 

For the non-factorizable annihilation diagrams Fig.(lg) and (Ih), all three wave functions 
are involved in the amplitudes. The integration of 63 can be performed by the 6 function 
5(62 — 63). The expressions of contributions of these two diagrams are 

32 f^ [^ 

M-aT = — CpT^rn^B / dxidx2dx3 I bidbib2db2(j)B{xi,bi) 
3 Jo Jo 

X { [-r-Tr-o {<I)p{x2){<I)t{^3){x2 - 1 + 0:3) - 0r(x3)(a;2 - 1 - X3)) 

+0p(a;2) (0^(3^3) (1 -X2- X3) + 0r(^3)(a;2 - X3 + 3))) + 0p(x2)0T(a;3)a;2] 

■hanfl{Xi,X2, X3, 61, b2)Eanf{tg) 

+ [r-Tro (0p(x2) (0^(2^3) (3^2 - X3 + 1) + 0^(2^3) (3^2 + 3:3 - 1)) 
-0p(a:2)(0T(2^3)(a;2 - 3:3 + 1) + 0^(3:3) (0:2 + X3 - 1))) 

+0p(a;2)0T(a;3)(a:3 - 1)] hanf2{xi, X2,X3, bi,b2)Eanf{th)} , (26) 



32 f f°° 

A4aT — 'irC! FT^rn^B / dxidx2dxs / bidbib2db2(f)Bixi,bi) 
3 Jo Jo 

X { [rT(j)i{x2){(j)T{x3) - (f)Tix3)){x3 + 1) - ro^Tl^a) (^^(xz) + 0p(a;2)) 

■ (X2 - 2)] hanfl{Xi, X2, X3, 61, b2)Eanf{tg) 

+ [ro0T(a;3)a;2(0p(a;2) + 0?(a;2)) - rT<pj>{x2){(t)T{.X:i) - 0r(a;3))(x3 - 1)] 

• hanf2{Xi, X2, X3, 61, b2)Eanf{th)] , (27) 

32 [^ f°° 

J^aT = "^"C'irvrm^ / dxidx2dx3 / bidbib2db24>B{xi,bi) 
3 Jo io 

X { [-rTro0p(x2)(0T(a;3)(a;2 - 1 + X3) + 4>t{x3){x2 - 1 - 0:3)) 
+rorT0p(x2) (07(3^3) (3^2 - X3 + 3) + 4>t{x3){x2 + X3 - 1)) 

+0^(x2)0T(a;3)(a;3 - 1)] hanfl{Xl,X2,X3,biM)Eanf{tg) 

+ [-^O^T0p(a;2)(0T(^3)(^2 + 1 - X3) + 0^(X3)(1 - X2 - X3)) 

-rorT0|^(x2)(0Ua;3)(-a;2 + 0:3 - 1) + 0^(2^3) (2^2 + 2:3 - 1)) 

+0p(x2)0r(x3)a;2] hanf2{Xl,X2,X3MM)Eanf{th)] ■ (28) 

If we exchange the pseudoscalar meson and the tensor meson in Fig.l, the result will be 
different. Because a tensor meson can not be produced through (V ± A) or tensor current, 
the factorizable emission diagrams do not contribute to the amplitude of B decays with a 
tensor meson emitted [3", '4l]. Therefore, there are only six diagrams left shown in Fig. 2. The 
individual decay amplitudes for these diagrams can be easily deduced from eg. ( 12U1I2S]) by 
the replacement of the wave functions of the pseudoscalar and the tensor meson, 

(t)j,{x) -^ -<t)T{x), (t)l{x) -^ </.^(x), <fp{x) -> 0^(x), 

0r(x) ^ -0^(x), 0^(x) ^ -0^(x), 0^(x) ^ -<fp{x), (29) 

tt -^ ro, To -^ tt. 

In addition, we must add a minus sign to M^rf after applying the above replacement. 

For the 39 -B — )■ PT decay channels, not all the effective operators contribute to each 
decay mode. We list the number of effective operators contributing to the individual decay 
channels in the Appendix B for reference. 
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FIG. 2: Diagrams contributing to the B — )• PT decays, with a tensor meson emitted. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

For the numerical analysis, we need various input parameters, such as decay constants, 
CKM elements and the wave functions, which are given in Appendix A. The CP-averaged 
branching ratios for those B — )■ PT decays with AS* = 1, together with Isgur-Scora- 
Grinstein-Wise II (ISGW2) model 2J] and the QCDF results [^ are shown in table I. 
The experimental data are taken from Ref.[l| and Ref. 37|. Similarly, the branching ratios 
of -B — 7- PT decays with AS* = calculated in the PQCD approach are shown in Table II. 
For illustration, we classify these decays to categories by their dominant topologies indicated 
through the symbols T (color-allowed tree), C (color-suppressed tree), P (penguin emission) 
and PA (penguin annihilation). Although we include also the W annihilation and W ex- 
change diagram contributions, none of these channels has dominant contribution from these 
two topology. For the theoretical uncertainties in our calculation, we estimated three kinds 
of them: The first errors are caused by the uncertainties of the decay constants of tensor 



11 



mesons; The second errors are from the decay constant {Jb = ( 0.21 ± 0.02) GeV)of B meson 



28l 



m 



and the shape parameter (ub = (0.5 ± 0.05) GeV)in the B meson wave function 
The third errors are estimated from the unknown next-to-leading order QCD corrections and 
the power corrections, characterized by the choice of the Xqcd = (0.25 ± 0.05) GeV and 
the variations of the factorization scales shown in Appendix B, respectively. It is easy to see 
that the dominant errors for the branching ratio calculations are from the non-perturbative 
wave functions. 

It is easy to see that there are large theoretical uncertainties in any of the individual decay 
mode calculations. However, we can reduce the uncertainties by ratios of decay channels. 
For example, simple relations among some decay channels are derived in the limit of SU(3) 
flavor symmetry 



1 



*o^+\ 



B{B+ -^ k;% 

i3(50 -^ a^K+) B{B+ -^ a+K^ 



(30) 



B{B+ -^ aOir+) B{BO -^ a^iTO) 

One can find that our results basically agree with the relation given above within the errors. 
Among considered B — t- PT decays, the PQCD predictions for the CP-averaged branching 
ratios vary in the range of 10^^ to 10~®. From the numerical results, we can see that the 
predicted branching ratios of penguin-dominated B — t- PT decays in PQCD are larger than 
those of naive factorization 2J, |25|, |27|] by one or two orders of magnitude, but are close to 
the QCDF predictions [^. For the leading tree-dominated modes such as a^vr"*" and /2vr"'", 
the predicted results in PQCD are bigger than those obtained by QCDF [4] but smaller 



than Ref. 22 ■ The reason is that the B to tensor form factor in this work is larger than that 



used in Ref. J]. But for a^^"*", the result is not larger than but the same as Ref. [4^]. This is 
the result of destructive interference from other topologies. It is worth of remarking that 
B^ — 7- K^^K^ and B^ — t- K2~K~^ are pure annihilation modes, which can be perturbatively 
calculated in the PQCD approach. 

The decays with a tensor meson emitted are prohibited in the naive factorization approach 
for the reason that a tensor meson can not be produced from the local (V ± A) or tensor 
currents 3|,|j]. In order to predict these decay channels, it is necessary to go beyond the naive 
factorization framework to estimate the contributions of the nonfactorizable and annihilation 
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TABLE I: The PQCD predictions of CP-averaged branching ratios (in units of 10 ^) for B — )■ PT 
decays with AS = 1, together with Isgur-Scora-Grinstein-Wise II (ISGW2) model [24] and QCDF 



results [J]. The experimental data are from Ref.[l[ and Ref. 



:Q 



37|. 



Decay Modes class 



This Work 



ISGW2 [24] QCDF [4] 



Expt. 



B+ -^ Kf 7r+ 


PA 


0.9 


B+ -^ K;+7r^ 


PA 


0.4 


B+ -^ 4k+ 


T,PA 


2.1 


B+ ^ a+i^o 


PA 


3.1 


B+ -^ f2K+ 


T,PA,P 


11.8 


B+ -^ f'K+ 


P,PA 


3.8 


B+ -^ k;+7] 


PA,P 


0.8 


B+ -^ K;+7]' 


PA,P 


12.7 


B^ -^ K*+7r- 


PA 


1.0 


50 ^ ^*o^o 


PA 


0.6 


B° -^ a^K+ 


T,PA 


5.0 


B^ ^ alK^ 


PA 


2.0 


50 ^ /2K0 


PA,P 


9.2 


B^ -^ tlK"^ 


P,PA 


3.7 


B^ ^ KfT] 


PA,P 


1.0 


B^ -^ Kfrj' 


PA,P 


11.6 



+0.2 +0.2 +0.3 




-0.2-0.2-0.2 




+0.1+0.1+0.1 
-0.0-0.1 -0.1 


0.090 


+0.7+0.6+0.6 
-0.6-0.5-0.5 


0.31 


+0.9+0.9+1.1 
-0.8 -0.8 -0.9 


0.011 


+2.7+3.2+3.0 
-2.4-2.8-2.7 


0.34 


+0.4+0.9+1.0 
-0.4 -0.8 -0.8 


0.004 


+0.2+0.3+0.3 
-0.2-0.2-0.3 


0.031 


+3.7 +4.5 +4.0 


1.41 


-3.2 -3.5-3.5 


+0.2 +0.2 +0.3 




-0.2 -0.2 -0.2 




+0.2+0.1+0.2 
-0.1 -0.1 -0.1 


0.084 


+1.6+1.4+1.3 
-1.4-1.1 -1.0 


0.58 


+0.5 +0.4 +0.6 
-0.5-0.4-0.5 


0.005 


+2.0+2.5+2.6 
-1.8-2.1 -2.2 


0.005 


+0.3+0.7+0.9 
-0.4-0.8-0.9 


0.00007 


+0.2+0.3+0.3 
-0.2-0.2-0.3 


0.029 


+3.6 +4.2 +3.8 
-2.9 -3.1 -3.1 


1.30 



3.1: 

2.2 



h8.3 

-3.1 

+4.7 

-1.9 



12.4 



+21.3 
-12.4 



5.6 



+2.2 
-1.4 



4 g+8.4 


<45 


8At]'^ 




q + 7.8 

-^-O-S.O 


1 06+0-28 
l.UD_Q29 


^■otli 


< 7.7 


6.81^3^^ 


9.1 ±3.0 


•-^■-^-12.1 


28.0l|3 


q q+8.5 
0.0_32 


<6.3 


-1 0+4.3 
-^■^-1.3 


<4.0 


^•'-8.1 




A 0+8.3 

^■^-3.5 




3Atli 


9 7+1.3 
^•'-1.2 


3.810 




6.6+if 


9.6 ±2.1 



13.7 



+3.2 
-3.1 



diagrams. Fortunately, in the PQCD approach, the contributions of the nonfactorizable 
diagrams with a tensor meson emitted (Fig. 2, c and d) are sizable and larger than that of 
the nonfactorizable diagrams emitting a pseudoscalar meson (Fig.l, c and d). The reason 
is that the asymmetry of the light-cone distribution amplitudes of the tensor meson makes 
the contributions from Fig. 2(c) and (d) strengthen with each other, while the situation is 
contrary for Fig. 1(c) and (d). One can see from Table II that for B — )■ 02 tt decays, the 
a^TT^ and a^vr" modes are highly suppressed relative to a^7r+ and 02^+, respectively. It is a 
natural consequence of factorization as the tensor meson can not be created from the (V-A) 
current. For B — t- a27r+(a2"7r+), the dominant contribution is from color-allowed factorizable 
emission diagrams, while for B — t- a^7r°(a^7r~), this large contribution is prohibited for the 
above reason. Therefore for 5+ — )■ a^Tr^, the left factorizable emission diagrams are color- 
suppressed, and for B^ — t- a^vr^, the dominant contribution is from nonfactorizable emission 
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TABLE II: The PQCD predictions of CP-averaged branching ratios (in units of 10^^) for B — t- PT 
decays with AS = 0, together with Isgur-Scora-Grinstein-Wise II (ISGW2) model |24;] and QCDF 
results [4|. The experimental data are from Ref.[l| and Ref.[37|. 



Decay Modes 


class 


This Work 


ISGW2 [24] 


QCDF [4] 


Expt. 


B+ -^ a07r+ 


T 


r,Q -I +12.8+14.2+3.1 
^y-J^ -10.6-10.4-2.8 


26.02 


301- 




B+ -^ a+vrO 


T,C 


n q +0.0 +0.1 +0.0 
"-'•"-' -0.0 -0.1 -0.0 


0.01 


2Atii 




B+ -^ a^r] 


C,PA,P 


1 n +0.3 +0.4 +0.4 
■L-U _o.3 _o.3 -0.3 


2.94 


1 1+2-8 
-■-•-L-l.l 




B+ -^ a+i 


C,PA,P 


q r +1.4+1.6+1.1 

•J-^ -1.0 -1.1 -0.8 


13.1 


1 1+"^-'^ 
-'-•^-1.2 




B+ ^ 


/2vr+ 


T 


An r +18.9 +18.9 +4.2 
^^•^ -15.4-13.9-3.9 


28.74 


r,7+14 


15.7li^ 


B+ ^ 


/^vr+ 


T 


1 9 +0.3+0.4+0.1 
-■-■^ -0.2 -0.3-0.1 


0.37 


omtrot 




B+ ^ 


k;+k^ 


PA,P 


1 9 +0.2 +0.2 +0.3 
^•^ -0.2 -0.2 -0.3 


4.0 X 10"^ 


4Atli 




B+ ^ 


KfK+ 


PA 


n o +0.1+0.2+0.3 
'-'•O -0.1 -0.2-0.2 




1 2+5-2 




B° -^ a^7r+ 


T 


Qo Q +35.1 +42.6 +5.8 
yo.» _29 9 _32.o -9.7 


48.82 


co+lS 


<3000 


B° -^ a+7T- 


T,PA 


9 7+0.5 +0.8+0.4 
^•' -0.3-0.5-0.3 




9 1+4.3 




B^ -^ aOyrO 


C 


^ f; +1.2 +1.6 +0.9 
^•" -1.0 -1.2 -0.7 


0.003 


2 4+^^-2 




fiO ^ air] 


C,PA,P 


r, f. +0.1 +0.2 +0.1 
"-*•" -0.1 -0.1 -0.1 


1.38 


6+1-6 




5° ^ ali 


C,PA,P 


1 o +0.6+0.7+0.4 
^•° -0.5 -0.6-0.4 


6.15 


n 5+2-2 

U.0_o.4 




fiO^ 


/2^° 


C 


9 o +0.7+0.7+0.6 
^•° -0.6 -0.6 -0.4 


0.003 


1 ^+4-2 




fiO^ 


/^vrO 


P 


n 9 +0.0 +0.1 +0.0 
'-'•^ -0.0-0.1 -0.0 


4.0 X 10"^ 


^•'-'^-0.05 




5°^ 


f2V 


C,P,PA 


9 f- +0.7 +0.8 +0.7 
^•0 -0.5 -0.6 -0.6 


1.52 


1 7+2.3 
-*■• ' -1.2 




B^^ 


f2V' 


C,PA,P 


q q +1.0+1.1 +0.9 
•^••^ -0.8 -0.9 -0.9 


6.8 


1 0+2:2 
-■^•"^-1.3 




B^^ 


f^V 


PA,P 


n nsj +0.03 +0.03 +0.01 

^■^° -0.02 -0.03 -0.02 


0.02 


0.021°:°^ 




fiO^ 


fW 


PA,P 


n no +0-00 +0-02 +0.02 

"-"•"-"^ -0.00 -0.02 -0.03 


0.09 


omt^l 




fiO^ 


k;+k- 


PA 


n ic +0.02+0.03+0.03 
•J--"^" -0.03-0.04-0.03 




n 3+0-'^ 

U.c5_o,2 




fiO^ 


k;-k+ 


PA 


r, Q +0.1 +0.3 +0.2 
^■^ -0.1 -0.1 -0.2 




1 0+1.6 

-■■•"^-1.0 




fiO^ 


KfK^ 


P,PA 


1 c +0.3 +0.3 +0.5 
^•^ -0.3-0.3-0.4 


3.0 X 10-4 


c 4+8.8 




B^-^ 


KfK^ 


P,PA 


n 0+0.1 +0.2+0.3 
'-'•0 -0.1 -0.1 -0.2 




9 9+5.4 

^ .^ 9 9 





diagrams suppressed by Wilson coefficient Ci. 

From table IVIt one can see that the factorizable contributions for the 5+ — ?■ 7^2°^^ 
and B^ -^ i^2^7r- decays are because of the emitted meson in these diagrams is the 
tensor meson. The contributions from nonfactorizable diagrams are suppressed by the small 
Wilson coefficients C3 and C5. Therefore the dominant contribution comes from the penguin 
annihilation diagrams. From table HI one can see that our predictions for the 5+ — > i^2°^^ 
and B^ — )■ i^2^7r- decays are much smaller than that from Ref. [^. The reason is that in 
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Ref.j^, there is an extremely large contribution from the quark loop diagrams. In PQCD 
approach, the quark loop correction is next-to-leading order and not considered in this 
work. In the B — )■ /2-ft' decays, we have tree diagram contribution as well as penguin 
emission diagram contributions, thus makes the branching ratios much larger than that of 
the B^ — ^ K2^'K^ and B^ -^ K'^'^-k^ decays. The current experimental measurements still 
have very large error bars. We expect the future experiment to give more information for 
these decays. 

For B -^ ir^r/W and B -^ as?]^') decays, one finds that B{B -^ K^r]') > B{B -> K^-q) 
and B{B — )• 02?]) ^ B{B — )■ a2rj'). For these modes, both rjq and rjg will contribute, but the 
relative sign of the rjg state with respect to the rjq is negative for the 77 and positive for the 
77', which leads to a destructive interference between rjq and r/^ for B — )• K2T1 and B — )• 02?], 
but a constructive interference for B — )■ K2T1' and B — )■ 02?]'. This is very similar to the 



situation for B -^ Kr]^'^ and B^ -^ K+tj^'^ decays [Ssl, l39|. 



TABLE III: The PQCD predictions of direct CP asymmetries(%) for B ^ PT decays with AS" = 1, 
comparison with the QCDF results J]. The experimental data are from Ref.[l|]. 



Decay Modes 



This Work 



QCDF [4] 



Expt. 



B+ -^ KfTT+ 


c t; +0.3 +2.6 +1.6 
0-0 -0.4-0.0-1.2 


1 fi+2-2 


B+ -^ k;+tt^ 


p. Q +2.6 +1.6 +3.7 
D.y _2,9_i,i _3.6 


n 9+17.8 
U-^-14.8 


B+ -, 4k+ 


c-r, Q +2.0 +2.1 +8.6 
J^-» _2.2 -0.4-10.1 


07 1 +33.3 

^'•-^-35.0 


B+ ^ a+K^ 


9 Q +0.1 +0.1 +0.5 
^•^ -0.1 -0.2 -0.8 


-o.6l°:| 


B+ -^ f2K+ 


94 fi +1-5 +2.4 +4.6 
^^•D -1.0 -2.6 -5.9 


-39.5+^1^ 


B+ -^ f'K+ 


c +1.5 +1.4 +1.5 
"•O -1.6 -1.0-1.8 


-n fi+^-3 


B+ -^ K;+r] 


cr /. +1.1 +2.2 +2.3 
^•■^ -0.6 -2.0 -1.3 


-^•^-5.6 


B+ -^ K;+r,' 


9 +0.1 +0.1 +0.9 
^•^-0.1 -0.3-0.5 


-, 7+3.2 
--■-•' -3.9 


B^ -^ k;+tt~ 


17 c: +1.4+1.6+2.7 
-^'•^ -1.6 -1.8 -1.3 


1 7+^-2 


BO ^ Kfn^ 


-in 7+0.1+1.7+1.9 
'''^■' -0.0-1.8-1.8 


7 1+23.5 
'•^-24.1 


B^ -^ a^K^ 


AQ q +1.9+1.3+7.1 
'iO.o _2 4 _o.3_9.9 


91 c;+28.9 


B^ ^ 4k^ 


-, Q +0.5 +0.4 +0.6 
-■-•y -0.5 -0.4 -0.5 


fi 7+6.5 


B'> ^ /2i^o 


-I 9 +0.3 +0.5 +0.2 
-L-^ -0.2-0.5-0.1 


-7.31^,:^ 


5° ^ /^i^o 


1 r, +0.1 +0.0 +0.0 
^■^ -0.3-0.1 -0.1 


o.s^J:? 


B^ -^ Kfr] 


c- n +0.5 +0.2 +1.7 
3-^-0.4 -0.1 -1.7 


Q 9+16.5 

-^•^-4.8 


50 ^ j^*0^> 


n. 7+0.1+0.1 +0.3 
"•' -0.0 -0.0-0.2 


9 9+3.3 
~^-^-4.0 



-+29 
'-24 



-es.oli? 



-45 ± 30 



-7.0 ±19.0 



We also give the direct CP asymmetry parameters for those B — t- PT decays with AS = 1, 
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together with the QCDF results ^ shown in table IIIII The experimental data are taken 
from Ref.pJ]. Similarly, the direct CP asymmetry parameters of i? — t- PT decays with 
AS = calculated in the PQCD approach are shown in Table IIVI The origin of theoretical 
uncertainties shown in these two tables are the same as those of the branching ratios in 
table H] and [ITl However, the dominant uncertainty here is the third one from the unknown 
higher order QCD corrections, since the hadronic parameter uncertainty mostly cancels due 
to the fact that the CP asymmetry is defined as the ratio of branching ratios. 

It is easy to see that some channels have very large direct CP asymmetries. But many of 
them have small branching ratios to make them difficult for experiments. We recommend the 
experimenters to search for the direct CP asymmetry in the channels like B^ — )■ j^K^, B^ — )■ 
a2K^, B^ —7- a2i]' and B^ — t- /27r+, for they have both large branching ratios and direct 
CP asymmetry parameters. In fact, there are already some experimental measurements 
for the CP asymmetries shown in table IIIII and IIVI Although the error bars are still large, 
we are happy to see that all these measured entries have the same sign as our theoretical 
calculations. This may imply that our approach gives the dominant strong phase in these 
channels. The decays 5°(5°) -^ a^n+fa+n-, B^{B^) -^ K;+K-/K;~K+ and B^{B^) -^ 
X*^K^/K2^K^ have a very complicated CP pattern through the B^B^ mixing. Four decay 
amplitudes are involved for each group of decays with 5 CP parameters to measure. We 



refer the readers to the similar situation for B^{B^) — > p 7r+/p+7r decays 40 1. 

For the decays involving fo in the final states, we have taken the f2~f2 mixing (Eq. (IA13|) ) 



.a 



into account, while in Ref. j^, /2 is considered as an {uu + dd)/y/2 state and /g a pure ss 
state. Although the mixing angle is small, the interference between /| and /| can bring 
some remarkable change. For example, the branching ratio of B^ — )■ /27r+ is bigger than 
the prediction in Ref. j^. This can be understood as follows: Because of the contribution 
from the color-allowed factorizable emission diagrams, although suppressed by the mixing 
angle, the contribution of /| term is at the same level with that of /| term. Due to the 
enhancement from /| term, the branching ratio becomes larger than the prediction without 
taking the mixing into account. The mixing can also bring remarkable change to direct CP 



B 



asymmetry. For B — t- f2v'''\ the direct CP asymmetries are zero jj] when /2 is a pure ss 
state. Since the direct CP asymmetry is proportional to the interference between the tree 
and penguin contributions [30], it should be zero indeed because there are no contributions 
of penguin operators when /g is a pure ss state. When taking the mixing into account, /| 
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TABLE IV: The PQCD predictions of direct CP asymmetries(%) for B ^ PT decays with A5 = 0, 
comparison with the QCDF results |4]- The experimental data are from Ref.[l|]. 



Decay Modes 


This Work 


QCDF [4] 


Expt. 


B+ ^ 


4tt+ 


n fi +0-1 +0.4+0.2 
"•0 -0.1 -0.5 -0.6 


9.61^^1 




B+ ^ 


a+^0 


r o +0.1+21.3+75.8 
""••O -0.1 -12.4-44.7 


-24.31^^,^ 




B+ ^ 


a+7? 


on Q +8-4 +9.6 +12.3 

-yu.y _3 7_^Q_5^ 


97 fi+'^3.4 
^' -0-127.6 




B+ ^ 


a+T]' 


A A K +0.8 +1.3 +6.8 
^^•^ -0.5 -0.2 -8.8 


q-| q+61.3 
"-•-^•"^-131.3 




B+ ^ 


/2vr+ 


97 R +3-4 +1-0 +8.9 
^'•" -2.5-1.4-7.1 


60.2trli 


41 ±30 


B+ ^ 


/^vr+ 


n ni. +0-1 +9-6 +13.8 
'^■'^^ -0.1 -8.9 -15.8 


0.0 




B+ ^ 


k;+k° 


An 7+1.3+1.8+16.4 
^'^•' -2.0-0.5-12.4 


30.3_33;7 




B+ ^ 


KfK+ 


4Q c; +4-7 +3.1 +23.5 
'iif.o _4 2 _4.8 -13.1 


-0.26l°:f, 




B^^ 


aOvfO 


c:q c; +4.7 +6.9 +4.2 
^•^•O -3.8-6.9-3.5 


-86.21128^9 




B^^ 


a^r] 


-17 7+17.7+11.2+21.8 
^'•' -15.7-22.6-24.5 


7fi 7+10° 




B^^ 


a\ri 


CQ Q +0.6 +10.0 +7.2 
^-'•y -0.0-6.0-7.0 


-66.0tM\ 




B^^ 


/2vr0 


Q o +13.9+2.8+11.8 
^•° -13.2 -7.5 -10.8 


07 n+103.8 
~'3'-^-85.5 




B^^ 


/^vrO 


n 7+2.7+1.0+6.8 
"•' -2.5-1.8 -6.4 


0.0 




B'>^ 


h-n 


49 rr +1.7+1.4+9.1 
'iz.o _-^ -j^ _i.8_9,8 


fiq 7+25-7 
"y-' -102.7 




B^^ 


hrl 


n nt^ +0-1 +5-0 +5.3 
U-U^ -0.6-5.1 -5.3 


09 q+22.9 

o^-<J-94.8 




B^^ 


f^r] 


7n Q +0.0+11.0+11.0 

m.y -2.7-15.2-12.3 


0.0 




B^^ 


fW 


4c; c; +3.2 +13.5 +18.5 
^'^•^ -6.8-12.1 -18.8 


0.0 





term can provide penguin contributions, then the direct CP asymmetries are no longer zero 
in this work. 

For S — > /2 1]^'' and /g r]^'' decays, the relevant final state mesons contain the same 
components -j={uu + dd) and ss, therefore they have the similar branching ratios. The small 
differences among their branching ratios mainly come from the different mixing coefficients, 
i.e., COS0, sin0, cos 6* and sin 6' (see Appendix A). 

V. SUMMARY 

We studied the charmless hadronic B — t- PT decays by employing the PQCD approach 
based on the kj- factorization. In addition to usual factorization contributions, we also calcu- 
lated the non-factorizable and annihilation type diagrams. From our numerical calculation 
and phenomenological analysis, we found the following results: 
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• The factorizable amplitude with a tensor meson emitted vanishes because a tensor 
meson cannot be created from the {V ± A) or {S ± P) currents. For these decay 
modes, the non- factorizable and annihilation diagrams' contributions are important. 
For example, B^ — ?► i^2'^7r+ and B^ — )■ ii'2^7r^ have sizable branching ratios because 
of the contributions of penguin annihilation diagrams. 

• For penguin-dominated B -^ PT decays, because of the dynamical penguin enhance- 
ment, the predicated branching ratios are larger by one or two orders of magnitude 
than those predicted in the naive factorization approach but close to the QCD factor- 
ization predictions in Ref. ^ 

• For tree-dominated decay modes, the branching ratios predicted by PQCD are usually 
very small except for a^^"*", a^vr"'" and /2vr"'" modes with branching ratios of order 10~^ 
or even larger. This basically agrees with the situation in Ref. j^ and Ref. [27 1. 

• For B — ^ K^T]^'^ decays, we find B{B — )■ -ft'a'?') ^ -^(-^ ~^ ^ll)- This large difference 
can be explained by the destructive and constructive interference between rjq and rjg. 

• From our calculation, we find that the interference between /| and /^ can bring some 
remarkable effects to some decays involving a /2 meson in branching ratio and direct 
CP asymmetry. 

• We predict large direct CP asymmetry for some of the B — )■ PT decays that accessible 
for the near future experiments. 
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Appendix A: Input Parameters And Distribution Amplitudes 

The masses and decay constants of tensor mesons are summarized in Table V. Other 
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TABLE V: The masses and decay constants of light tensor mesons 



Tensor(mass(MeV)) friMeV) /r(MeV) 



/2(1270) 
/2(1525) 
02(1320) 
K|(1430) 



102 ± 6 


117 ± 25 


126 ± 12 


65 ± 12 


107 ± 6 


105 ± 21 


118 ± 5 


77 ± 14 



input parameters are 



Ag = 0.25, m, = 4.8, U = 0.131, /^ = 0.16, 



ttIq = 1.4, mQ 



_ft: 



1.6, mg 



'7? 



1.07, mg^ 



1.92. 



(Al) 



We adopt the Wolfenstein parameterization for the CKM matrix, A = 0.808, A = 0.2253, 
p = 0.132 and 77 = 0.341 [Ij. 

The twist-2(3) pseudoscalar meson distribution ampUtude(s) 0p (0p, 0p) (P = vr, K) can 
be parameterized as 4l|, |42|] , 






2VQ 



2^6 
3/i^ 



x{l~x) 1 + 0.44C2/^(t) + 0.25Cf^(t) 

1 + o.43C2^/^(t) + omcy\t) 



2^6 
Ik 



x{l — x) 



2^6 L 



1 + 0.17Cf/^(t) + 0.2C2/^(t) 
1 + 0.24C2/^(t) - 0.11Cy\t) 
Cl^\t) + 0.35C3/'(t) 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 
(AT) 



The Gegenbauer polynomials can be defined by 



1 

2 
1 



■y3/2 



ct'\t) 



3t, 



■y3/2 



(3t^-l), C',^\t) = -(5t^-l), 



Cl'\t) 

C]/'(t) = i(35t^-30t' + 3), Cf'(t) = ^(21t^-14t2 + l), 



(A8) 



where t = 2x — 1. In the above distribution amplitudes for kaon, the momentum fraction 
X is carried by the "s" quark. 
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For the 77 — 77' system, we use the quark-flavor basis 43|], with rig and rjs deflned by 

rjg = —z={uu + dd), rjs = ss. (A9) 

v2 



The physical states 7/ and r]' can be given by 



cos — sin \ I rjq 
sin0 COS0 J \ Vs 



The decay constants are related to fg and fg via the same mixing matrix, 



fq fs 
J T) J r) 
-fQ fs 



COS (p — sm I 
sin cos d) 



U 
/. 



(AlO) 



(All) 



The three input parameters /„, fg and (j) have been extracted from related experiments 
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4^: 



fg = (1.07 ± 0.02)/,, /, = (1.34 ± 0.06)/,, = 39.3° ± 1.0°. 



(A12) 



Like the 77 — 77' mixing, the isoscalar tensor states /2(1270) and /2(1525) also have a 
mixing and can be given by 



/2 = /|cose + /^in^, 
/^ = /|sin^-/2^cos^. 



(A13) 



where /I = -js{uu + dd), /| = ss and the mixing angle 6 = 5.8° 45|, 7.8° 46| or (9 ± 1)° 

Q. 



Appendix B: Amplitude And Related Hard Functions 



For each individual decay channel, various effective operators contribute to the decay 
amplitude. We summarize the number of effective operators contributing to every channel 
in Table IVTl and NTH for the AS* = 1 and AS* = 0, respectively, with 



ai = — + 62, 02 = Ci + — , 
«j = Ci + ^(j = 3,5,7,9), a„ = ^ + C4n = 4,6,8,10). 



(Bl) 



20 



TABLE VI: The effective operators contributing 
cfianneis 



to each decay mode with AS* = 1 
annihilation 



emission 



factorizabie 



non-factorizabie 



factorizabie 



non- factorizabie 



50 - 

B+ 
B+ 

B+ 
B+ 
B+ 
B+ 
B+ 
B+ 



rfKf 
> Kf-K+ 



ai,a4)«6)fl8)«10 
04,06,08,010 

02,07,09 
04,06,08,010 

«3)«4) 05,00, 

OT^asjOgjOio 

04)06) ^SiOio 

02,07,09 

01,04,06,08)010 
01,04,06,08)010 
02)03)05)07,09 



03,04,05,06, 
07)08)09,010 



1)L'3, 05,07,09 

Cl)C'3)C'5,C'7)C'9 

2,03)05,07,08) 

C9,C'lO 

2,03,05,07,08) 

C9,Cio 
r^ r^ r^ r^ r^ 

02,03)04,05,06) 

C7,C8,C9,Cio 

C2,C'3)C'4,C'5)C'6) 

C7,C8,C9,Cio 
C3,C4,C5,C6,C7, 

C8,C9,Cio 

C3,C'4)C'5,C'6)C'7) 

C8,C'9)C'lO 

C3)C'5,C'7)C'9 

C3)C'5,C'7)C'9 

Cl,C'2)C'3,C'5)C'7) 

C8,C9,Cio 
Cl,C'2) (^3,(75,(77, 

Cs, (^9,(^10 
Cl,C'2)C'3,C'4)C'5) 
C6)C'7,C'8)C'9,C'lO 

1,02)03,04,05, 
r^ r^ r~i /^i r^ 

06)07,08) 09,010 

C3,C4,C5,C6, 
C7) (^8,(^9 ,(^10 
C3, (^4, (75, 6*6, 

C7)C'8,C'9)C'lO 



04,06 
04,06 

04,06 
04,06 
04,06 
04,06 
04,06 
04,06 



O8,Oi0 
O8,Oi0 

O8,Oi0 
O8,Oi0 
O8,Oi0 
O8,Oi0 
O8,Oi0 
O8,Oi0 



01,04,06)08)010 
01,04,06,08,010 
01,04,06)08)010 

01,04,06)08)010 

01,04,06)08)010 

01,04,06)08)010 

01,04,06)08)010 

01,04,06)08)010 



C3,C'5)C'7)C'9 
C3,C'5)C'7)C'9 

C3,C'5)C'7, C'g 

C3,C'5)C'7, C'g 
C3,C'5)C'7) C'g 

C3,C'5)C'7, C^ 
Cz.C^-.C'j, Cg 

C3,C5,Ct, Cg 
r^ r^ r^ r^ r^ 

01,03,05,07,09 

Cl)C'3,C'5)C'7,C'9 
r^ r^ r^ r^ r^ 

01,03,05,07,09 

ry ry /^ ry /^ 

01,03,05,07,09 

ry ry /^ ry r^ 

01,03,05,07,09 

r^ r^ r^ r^ r^ 

01,03,05,07,09 

ry ry /^ ry /^ 

01,03,05,07,09 

r^ r^ r^ r^ r^ 

01,03,05,07,09 



For factorizabie emission diagrams Fig.l. (la) and (lb), the h function is given by 



/ie/(a:i,X3,6i,63) 



X {6'(6i - 63)7^0 (ySs^^B^i) /o {\fximBbz) 
+9{bs - bi)Ko {^/x^mBh) Iq {y/x^mBh)} 
xSt{x3). 



(B2) 



The hard scales 



ta = max{y/x^mB, l/h, I/63} 
tb = max{y^mB, l/h, I/&3}, 



(B3) 



are the maximum energy scales in each diagrams to cancel the large logarithmic radiative 
corrections. The 5*^ re-sums the threshold logarithms In^ x in the hard kernels to all orders, 
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TABLE VII: The effective operators contributing to each decay mode with AS* = 

annihilation 



channels 



i?0- 
B^- 
B^- 
B^- 
50- 

50 - 

5+ 
5+ 

B+ 
B+ 
B+ 
B+ 






IT a, 



2 



fir]' 

fiv' 



a07r+ 



On vr 



/l^^ 



K+Kf 



emission 



f'actorizable 



02,04,06,07, 
02,03,04,05,06, 

07)^8) ^9)010 
01,04,06,08,010 



02,04,06,07,08, 
09,010 

a3,a5,«7,«9 
02,03,04,05,06, 

a7)a8,a9)Oio 
a3ja5)«7)«9 



O-i, 0-6,0-8 J O-IO 



01,04,06,08,010 

02,04,06,07,08, 
09,010 

01,04,06,08,010 

02,«3; 04, 05,06, 
07)^8)091010 

03)05)07,09 



04)06)08,010 



non- f'actorizable 

<-'2)'-'3)<-^4)<-^5)'-'6) 

C7)C'8)C'9)C'lO 
C2,C^,C4,C^,Cq, 

C7,C8,Cg,Cio 

r^ r^ r^ r^ r^ 

L'l)L'3)<-^5)L-7)<-^9 

r^ r^ r^ r^ r^ 

L'l)L'3)<-^5)L-7)<-^9 

C2)C'3,C5,C7,C8, 

C9)C'lO 

C4,C6,C8,Cio 

C4,C6,C8,Cio 

C2,C'3,C4,C5,C6, 
C7)C'8) (^9)^^10 

C4,C6,C8,Cio 
C4,C6,C8,Cio 



C3)C'5)C'7)C'9 

C3)C'5,C7,C9 

r^ r^ r^ r^ r^ 

<-^l)L-2)<-^3)<-^5)L'7) 

C8)C'9)C'lO 
Cl)C'2,C3,C5,C7, 

C8)C'9)C'lO 
Cl,C'2,C'3)C'4,C'5) 
C6)C'7)C'8)C'9)C'lO 

Cl,C'2)C'3)C'4)C'5, 
r~i /^ r~i r~i /^ 
L'6)<-^7)L^8)L^9)<-^10 

C4,C6,C8,Cio 

C4,C6,C8,Cio 

C3,C5,C7,Cg 



f'actorizable 



02,04,06,07, 

08)09,010 

02,04,06,07, 

08)09,010 
02,03,04,05,06, 

07)8 ) 09,010 
02,03,04,05,06, 

07,8 ,09,010 
02,03,04,05,06, 

07)8 ) 09,010 



02,03,04,05,06, 
07)08,09,010 
03)05)07,09 



02)03,05,07,09 
02,03,05,07,09 
03,04,05,06, 
07)08)09,010 
03,04,05,06, 
07)08)09,010 

01,04,06,08,010 
01,04,06,08,010 
01,04,06,08,010 
01,04,06,08,010 



01,04,06,08,010 
01,04,06,08,010 



non- f'actorizable 



C27C37C57C7! 

C8)C'9)C'lO 

C2)C'3)C'5)C'7) 

Cs,Cg,Cio 

C2,C'3,C4,C5,C6, 

C7)C'8)C'9)C'lO 
C2)C'3,C4,C5,C6, 

C7,C8,C9,Cio 

C2,C'3,C4,C5,C6, 

C7,C8,C9,Cio 



C2,C3,C4,C5,C6, 
C7,Cs,Cg,Cio 

C4 , Ce , 6*8 , Cio 



2,04,06,08,010 
C2,C4,C6,C8,Cio 
C3,C4,C5,C6, 
C7)C'8)C'9)C'lO 
C3)C'4)C'5)C'6) 
C7)C'8)C'9)C'lO 

r^ r^ r^ r^ r^ 

01,03,05,07,09 

Cl)C'3)C'5)C'7)C'9 

r^ r^ r^ r^ r^ 

01,03,05,07,09 

r^ /^ r^ r^ r^ 

01,03,05,07,09 



r~i f^ r~i r~i /^ 

01,03,05,07,09 

r^ r^ r^ r^ r^ 

01,03,05,07,09 



which is given by [33 1 



•S^tW = ^^,, , — ^[x(l-x)J 



(B4) 



v^r(i + c) 

with c = 0.3 in this work. In the nonfactorizable contributions, the St{x) provides a 
very small numerical effect to the amplitude j47j|. Therefore, we omit the St{x) in those 
contributions. 

The evolution factors Eef{ta) and Eef{tb) in the matrix elements (see section III) are 
given by 



Eef{t) = a,(t)exp[-^B(t)-53(t)]. 



(B5) 
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The Sudakov exponents are defined as 

Ssit) = s(xi-^,6ij + - — 7g(a,(/i)), (B6) 

S2it) = s (^2^,62) + s (^(1-X2)^,h^ + 2^*^ ^l,{asm, (B7) 

Ssit) = s (^3^,63) + s (^{l-xs)^,h^ + 2^*^ ^l,{asm, (B8) 



where the s{Q,b) can be found in the Appendix A in the Ref. 29 1 



For the other diagrams, the related functions are summarized as follows: 



tc = maxiy/xiXsUiB, a/|1 - Xi - X2\x3mB, I/61, 1/&2}, 

td = max{^/x^mB, a/|xi -X2|x3ms, I/61, 1/&2}, (B9) 

Eenf{t) = as{t) ■ exp[-SBit) - S^it) - Ssit)] I fe, = b3, (BIO) 

henf{xi,X2,X3,bi,b2) = [d{b2 - bi)Ko{y/x^mBb2)Io{y/x^mBbi) 

+9{bi - b2)Ko{y/xiX3mBbi)Io{y/^iX3'mBb2)] 

^-H'^q^ {J{x2 - Xi)x2,mBb2] , X2 - xi > 0; 

/ }- X ^ (Bll) 

Kq ( V (xi - X2)x3mBb2] , Xi - X2 > 0. 



te = max{Vl - X3mB, I/&2, 1/&3}, 

i/ = max{yi^mB, 1/62, l/fes}, (B12) 

Eaf{t) = as{t) ■ exp[-52(t) - Ss{t)], (B13) 



haf{x2,xs,b2,b3) = {—fH^^\y/x^mBb2) 



e{b2 - &3)^(S^^ {\^rnBb2) Jo (y/x^rriBbs) + 
^(63 - b2)H^'^ i^MniBbs) Jo (V^mB62)l ■ St{xs). (B14) 



tg = max{A/x2(l - X3)mB, a/1 - (1 - Xi - X2)mB, I/&1, 1/62} 

t/, = max{v^X2(l-X3)mB, a/|xi - X2|(1 - X3)mB, I/61, 1/62}, (B15) 

^an/ = as(t) ■ exp[-S'B(t) - 5'2(t) - S'3(t)] I 62=63, (B16) 
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hanfi{xi,X2,X3,bi,b2) = — 0{bi - &2)-f^o (^3:2(1 - X3)mBbi) Jo f A/a;2(l - X3)mBb: 

+6{b2 - bi)H^^^ i\/x^(l^^X3)mBb2) Jo iV^^i^ - X3)mBbij 
xKo f yr^^7r^^x7^^x^)a?^mB6ij , (B17) 

hanf2{Xl,X2,X3,bi,b2) = — 6{bi - ^2)^^^^^ (^2:2(1 - X3)mBbiJ Jo (^a/x2(1 - X3)mBh 

+6{b2 - bi)WQ f a/x2(1 -X3)mBb2J Jq f a/x2(1 - X3)mBbi 



f i^o ( a/(x2 -a;i)(l -X3)mBbi ) , xi - X2 < 0, 
X <; / ^ N ^ (B18) 

Ko [y {Xi - X2){1 ~ X3)mBbi\ , Xi-X2>0, 



where H^o\z) = M^) + ^^'ol-z)- 
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